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Abstract 

We develop a systematic method of the perturbative expansion around the 
Gaussian effective action based on the background field method. We show, 
by applying the method to the quantum mechanical anharmonic oscillator 
problem, that even the first non-trivial correction terms greatly improve the 
Gaussian approximation. 



*email: ghlee@ 'theory. yonsei. ac.kr 
^ email: jhyee@phya. yonsei. ac.kr 

1 



I. INTRODUCTION 



The variational Gaussian approximation method has provided a convenient and practical 
device in obtaining non-perturbative information from various quantum field theories, and 
has played an important role in enlarging our knowledge on the non-perturbative nature of 
the theories HI]. The main drawback of this method, however, is that it does not provide 
a systematic procedure to compute the correction terms to the approximation. One natu- 
ral way to improve the Gaussian approximation is to use more general non-Gaussian trial 
wave functional for the variational approximation B. Although this approach can easily 
be applied to the case of quantum mechanics, it has not been very successful in the case of 
quantum field theory. 

The efforts to establish the systematic approximation method based on the variational 
techniques have been made by authors of ref. and the so-called optimized expansion 
method and variational perturbation theory have been formulated. It is based on the fact 
that the Gaussian approximation is to approximate an interacting theory by a free field 
theory with a new mass, that is determined by variational gap equation. One constructs a 
complete set of energy eigenstates for this free theory and use it as a basis for perturbative 
expansion of the interacting theory. Cea and Tedesco used this variational perturbation 
theory to evaluate the second order correction terms to the Gaussian effective potential for 
the (1+1) and (2+1) dimensional 4> 4 theories. This perturbation theory is conceptually sim- 
ple and provides a very convenient systematic procedure to evaluate the correction terms to 
the Gaussian effective potentials for time-independent systems. For time- dependent systems 
such as those treated by Eboli, Jackiw and Pi Q, however, one needs a method that can 
deal with more general effective action. 

It is the purpose of this paper to develop a systematic method of computing correction 
terms to the Gaussian effective action based on the background field method(BFM) [|J. The 
background field method provides a very effective and simple way of computing effective 
action in the loop-expansion. The reason why the procedure of computing effective action 
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in BFM becomes simplified is that one has rearranged the summation of diagrams in such 
a way that the propagator used in the Feynman rule already contains the part of the 1-loop 
interaction effect. We proceed one step further so that the zeroth order term in the expansion 
would give the Gaussian approximation, and the propagator used in the expansion would 
be the one obtained in the Gaussian approximation. This provides a very convenient and 
simple method in computing order by order correction terms to the Gaussian approximation 
of the effective action. 

In the next section we give a brief review of BFM and develope the perturbative expansion 



method around the Gaussian effective action by using scalar (p A theory. In section [II], 
we illustrate the effectiveness of the method by evaluating the first non-trivial correction 
terms to the Gaussian effective potential of the (0+l)-dimensional (p 4 theory (anharmonic 
oscillator) and compare the result with those of various other methods. We conclude with 
some discussions in the last section. 

II. PERTURBATIVE EXPANSION AROUND THE GAUSSIAN EFFECTIVE 

ACTION 

For a system described by the classical action, 

S[<P] = Jd^xC^ix)^^)}, (1) 

the generating functional for Green's functions is defined by the path integral, 

< 0|0 > J = Z[J] = e lW[J] = J V$e iS ^+ l3 \ (2) 

where J(x) is the external source, W[J] is the generating functional for connected Green's 
functions, and the integral convention, J(f) = j d n xJ(x)(p(x) is used in the exponent. The 
Green's functions can be obtained by 

< O+ITM*!) • • • 0(s n )]|O_ > J = (— 6 —) ■ ■ ■ (^-^)e iW[J] . (3) 

ioJ(xi) ioJ(x n ) 
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The vacuum expectation value of the field operator in the presence of external source is 
defined by 

^ < *' > '=4 ,F|4 (4) 

The effective action is then defined by the Legendre transformation; 

r>] = w[j] - v , j. (5) 

Taking the derivative of T[ip] with respect to the external source J gives 

which is of the same form as the classical equation of motion. 

We now introduce a new action S[(f> + B] which is obtained by shifting the field by a 
background field B. This new action gives a new generating functional for Green's functions 
in the presence of the background field B, 

Z[J,B] = e lW ^ 



and the vacuum expectation value of field operator is now defined by 

0=jjW[J,B} = ^j\og(Z[J,B}). (8) 
Since Eq. (|7|) can be written as 

Z[J, B} = J vfeiSW+iJifi-B) = Z [J]e- iJB , (9) 

the new generating functional W[J, B] for connected Green's functions in the presence of 
the background field is related to the old one by, 

W[J, B] — W[J] — JB. (10) 

By taking the derivative of Eq. ( |T0"D with respect to J, one obtains the relation between </3 
and ip; 



<p = ip-B. (11) 

The new effective action T[<p, B] in the presence of the background field B can then be 
written as 

T[<p,B] = W[J,B] - Jip 

= (W[J) - JB) - J(ip - B) 

= T[<p + B]. (12) 

If we set if = 0, then Eq. (|12"D becomes 

T[B]=f[^ = 0,B], (13) 

which implies that the effective action T[B] consists of one c/5-particle irreducible diagrams 
with no external c/3-lines. Due to Eq. (|]) and (|12D, the effective action can be considered to 
be the extremum value of W[J, B] with respect to the variation of the external source J. 
We now consider the n-dimensional scalar 4 theory described by the Lagrangian density, 

L = -\^V-^-^\ (14) 

where V>~ x is defined by 

iV- 1 = (d^ + ^). (15) 

The generating functional for Green's functions in the presence of the background field 
B is given by 



Z[J, B) = J p0e i5( * +B)+ ^ J 



detV^Kexpl^BV^B - ^B 4 ] (16) 

x eMiBV -> _ exph « ( A-y _ p(^- 



where K is defined by 



K-^V^-^BXy. (17) 



The effective action is given by 

e ir ^ = Z[J,B]\^ =0 . (18) 

Note that the first line of Eq. fllTf ) gives rise to the 1-loop effective action and the last line, 
upon setting ip = 0, represents the higher loop contributions to the effective action. This 
higher loop contributions can be evaluated by using the Feynman rule with the propagator 
given by K xy of Eq. ([T7|). Thus the effective action consists of one (^-particle irreducible 
bubble diagrams with no external lines, and the procedure of evaluating the higher loop 
contributions is simplified compared to the conventional method. The reason why it is 
simplified is that the propagators used in the Feynman rule, Eq. flTTD, already contains the 
interaction effect through the background field B, and the generating functional has been 
rearranged in such a way that the zeroth order term in the expansion, the first line of Eq. 



fllq), is the one- loop effective action 

We now proceed one step further so that the propagator used for the perturbative ex- 
pansion would be the one obtained in the Gaussian approximation. In other words, we want 
to rearrange e* r ' B l in such a way that the zeroth order term in the expansion would give the 
Gaussian approximation of the effective action. To do this, we consider the relations: 

_L e y yGyzJz = J = (GJ) xe kJGJ (19) 

0J X 

1)2 e* J * G »' J ' = {{GJf x + G xx )e^ JGJ (20) 



^-e^ G ** J * = {{GJf x + 3G xx (GJ) x )e l * JGJ (21) 



5JI 
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,^ J ' <: ' ' = {{GJ)i + QG XX (GJ) 2 X + 3GL)e^ JGJ , (22) 

which would appear in the expansion of the last line of Eq. ([II]). Note that Eq. ([H|) and 
([22|) contain the diagrams where internal lines come out of a point and go back to the same 
point as shown in Fig. |1] and |2|, and the Gaussian approximation of the effective action 
consists of such diagrams. 

To extract such diagrams out of the last line of Eq. (|TED, we define new functional 
derivatives in such a way that 



63 "e l 2 JGJ = ((GJ) x ) 3 e k * JGJ (23) 



()i )' e * JGJ = ((GJ) x ) 4 e^ JGJ . (24) 



Comparing Eqs.([H]) and ( |22|) and Eqs. (|23|) and (p4]) , we find that the primed derivatives 
must be defined by 

"'' ~ 3G XX — — (25) 



sjr 6 J* 5 J, 



X 



<^S- 6G -4 +3GL (26) 

We note that the primed derivatives operated on eh JGJ do not generate the diagrams that 
contribute to the Gaussian approximation, such as those shown in Fig. 0(b), Fig.||(b) and 
(c). 

This implies that by using the primed derivatives one may extract the diagrams that 
contribute to the Gaussian effective action, out of the last line of Eq. fll6|). To do this we 
need to find a Green's function G which satisfies, 

eXPi 3! ( 5J j 4\ { 5J )le ~ 

Nexp i A h ]exp[ ir ( ^ }/ -¥ ( ^ ),]e|JGJ ' (27) 

where a normalization constant N and a function A(x) are to be determined. By using the 
definitions Eq. ( ff5] ) and (PB[), one easily find G, N and A that satisfy Eq. (fZTj): 

iV = , e^§GU, (28) 

A = ^BG B) (29) 
G X y = K X y + —G xx 5 xy 

= - jB 2 x S xy + jG x J xy . (30) 
Thus we can rewrite the generating functional for Green's functions, Eq. fll6|) as 
Z[J, B] = det v^Gexp^p-^ - ^S 4 + ^<£j 



x exp 



-iBV 1 — ^£? 3 + ^BG xx )jj 



exp[f(A) / -|(&) / ]e^, (31) 



where the functional derivatives in the last line do not generate the cactus type diagrams 
such as those shown in Fig.^(c). 

The effective action T[B] is given by 



e iT[B] = z[J,B]\ 



<p=0 



detV^eM^BV-iB-^ + jGU I[B], 



(32) 



where 



I[B] = exp 



XB 5 3 v iX. 5 A V1 i 



4! V <LL 4 ' 



JGJ 



(33) 



<p=0 



Note that Eq. (^), which determines the Green's function G xy , is exactly the variational 
gap equation written in Minkowski space notation IJ, and therefore the coefficient of I[B] in 
Eq. ( |32"D gives rise to the Gaussian effective action as we have alluded to above. Since I[B] 
of Eq. (p3|) can be expanded as a power series in A as is done in the conventional BFM, we 
have the perturbative expansion of the effective action around the Gaussian approximation. 

We now consider the contributions coming from I[B}. The linear term in the exponent 
of Eq. (|33| ) generate tadpole diagrams. This term does not contribute to the effective action 
because the effective action T[B] is the extremum value of W[J, B] as explained earlier, and 
the tadpole term does not change the extremum value. We therefore see that I[B] has the 
same structure as the last line of Eq. (|16|) except that the propagator is replaced by G xy , 
and the functional derivatives are replaced by the primed derivatives. Thus I[B] consists of 
one (^-particle irreducible bubble diagrams with no external c/3-lines, and without the cactus 
type diagrams (with internal lines coming out of a point and going back to the same point). 
We have thus established the systematic correction method to the Gaussian approximation, 
where one can compute order by order correction terms to the Gaussian effective action. 
One can easily show that the A 2 contribution of I[B] consists of the diagrams shown in 
Figfl 
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If we consider the space-time independent background field in the above formulation, we 
obtain the effective potential defined by 

V*M . -IM (34, 

Since B is a constant, G xy is a function of (x — y). Defining the Fourier transformation of 



G X y 



and the effective mass as 



G xy = I g(p)e«^ (35) 



m*=p? + ±&-$ lg(p), (36) 



we can write Eq. (|30|) as 



where 



2 2 J p 



2(jr — m z j 



(38) 



P , (27r)"' 

and n is the dimension of space-time. We note that Eqs.(^) and fl37|) are the well known 
mass gap equation for the Gaussian effective potential |l| . Thus the zeroth order contribution 
of Eq. (|32| ) is the Gaussian effective potential and the first non-trivial correction terms to 
the Gaussian approximation are of the order A 2 as shown in Fig.|3]. 

Thus up to the A 2 order contribution of I[B] the effective potential is given by 

V ef f = V G + V P , (39) 

where 

Va = ^ + A S 4 _ l_ {m 2 _ ^ _ A s2)2 + m lQg ^ (4Q) 

Vp ss - i^G\ (41) 

and G is given by Eqs.(^) ~ (ffTD- Vg is the Gaussian effective potential, and Vp is the first 
non-trivial contribution from I[B] which is shown in Fig.[3]. This result is the same as that 
of Cea and Tedesco [|J computed by using the variational perturbation theory. 



III. EFFECTIVE POTENTIAL FOR THE (0+1)-DIMENSIONAL 4 THEORY 



To illustrate how much our method improve the Gaussian approximation we now consider 
the (0+l)-dimensional 4 theory, which is the quantum mechanical anharmonic oscillator, 
and compare the results with those of various other methods. By applying the Wick rotation, 
Eqs.(|40|) and (f4l"D can be written in this case as 



where 



V eff [B] = V c + V P , (42) 



V g =£# + ±#-^-,?-\#r + df-, (43) 

V r =™# . (44) 

144m 4 1536 (y/irf) 5 

The effective mass, m, in Eqs. (|43"|) and fll!]) is given by 

m 2 = ^ 2 + ±B 2 + -^=, (45) 

where we have put the Plank constant explicitly. 

The improved Gaussian effective potential contains all diagrams up to 2-loop. Therefore, 
The 2-loop effective potential, which can be obtained by taking the terms up to h 2 order 
from Eq. (§33), is 



+ h 2 ( -(fi 2 + -B 2 )- 1 - — S 2 (/i 2 + -B 2 )- 2 ] . (46) 
V32 VP 2 ; 144 VP 2 ; J K ' 

Ignoring the ft 2 term, we obtain the 1-loop effective potential, which agrees with that of 

Ref . . Using these results, we evaluate the ground state energies of various approximation 

methods and compare their results. 

For simplicity, we consider the case of positive fi 2 — 1, and compare the ground state 

energies of various approximation methods in the unit h = 1. Then the classical potential 

has minimum V c i = at B = 0. The values of the coupling constant will be chosen to be 

that of 0. (A of H corresponds to our 4.) 
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Table | shows the ground state energy values of various approximations and their errors 
compared to the numerical results, for several coupling constants. This shows that the first 
non-trivial correction to the Gaussian approximation greatly improves the result for the 
small coupling region. For strong coupling, Gaussian and improved Gaussian approximation 
display much better results than the loop expansion, which shows the non-perturbative 
nature of the approximation methods. Though the GEP is accurate to within about 2% in 
the strong coupling region, the error of the improved GEP is smaller than about 0.8%. This 
shows that the perturbative expansion method greatly improves the Gaussian approximation 
even at the first non-trivial correction level. 



IV. DISCUSSION 

We have developed the perturbative expansion method around the Gaussian effective ac- 
tion by using the background field method. The method is based on the observation that the 
Gaussian effective action consists of cactus type diagrams, which were extracted out of the 
functional derivative part of the effective action, i.e., the last line of Eq. (|i6|), by introducing 
the primed derivatives defined in Eqs.(|25|) and (p6|) . This procedure effectively rearranged 



the diagrams in such a way that the propagator used in the perturbative expansion becomes 
that of Gaussian approximation, and the zeroth order term of the effective action consists 
of the Gaussian one. 

In the last section we have considered the quantum mechanical anharmonic oscillator, and 
have shown that the perturbative correction greatly improves the Gaussian approximation 
even at the first non-trivial (A 2 contribution) correction level. 

One can easily compute the A 2 contributions to the effective potentials for the (1+1) and 
(2+1) dimensional (p A theories, and show that they agree with those of Cea and Tedesco 
0. Our method appears to be simpler than the variational perturbation theory in practical 
applications and can also be used in the cases of time- dependent systems. Our method can 
easily be generalized to the cases of fermionic and gauge field theories. 
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Figure 1. Feynman diagrams for Eq. (f2~l|). 
Figure 2. Feynman diagrams for Eq. (^). 

Figure 3. Digrams contributing to the A 2 -order corrections to the Gaussian effective action. 
Table 1. ground state energy for anharmonic oscillator. 
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TABLES 

TABLE I. ground state energy for anharmonic oscillator 



A 


1-loop 


2-loop 


Gaussian 


improved Gaussian 


Eq 


error (%) 


E 


error(%) 


Eq 


error(%) 


Eq 


error (%) 





0.5 





0.5 





0.5 





0.5 





0.24 


0.5 


-1.43 


0.508 


0.048 


0.507 


0.006 


0.507 


-0.0003 


2.4 


0.5 


-10.6 


0.575 


2.84 


0.560 


0.21 


0.559 


-0.037 


24 


0.5 


-37.8 


1.25 


55.5 


0.813 


1.09 


0.801 


-0.37 


240 


0.5 


-66.8 


8. 


432 


1.53 


1.75 


1.49 


-0.68 


2400 


0.5 


-84.0 


75.5 


2311 


3.19 


1.95 


3.11 


-0.79 


oo 










0.236A3 


2.01 


0.230X3 


-0.82 
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